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ABSTRACT 

We present observations of the Zeeman effect in OH thermal absorption main 
lines at 1665 and 1667 MHz taken with the Karl G. Jansky Very Large Array 
(VLA) toward the star forming region S88B. The OH absorption profiles toward 
this source are complicated, and contain several blended components toward a 
number of positions. Almost all of the OH absorbing gas is located in the eastern 
parts of S88B, toward the compact continuum source S88B-2 and the eastern 
parts of the extended continuum source S88B-1. The ratio of 1665/1667 MHz OH 
line intensities indicates the gas is likely highly clumped, in agreement with other 
molecular emission line observations in the literature. S88-B appears to present 
a similar geometry to the well-known star forming region Ml 7, in that there 
is an edge-on eastward progression from ionized to molecular gas. The detected 
magnetic fields appear to mirror this eastward transition; we detected line-of-sight 
magnetic fields ranging from 90-400 /iG, with the lowest values of the field to the 
southwest of the S88B-1 continuum peak, and the highest values to its northeast. 
We used the detected fields to assess the importance of the magnetic field in S88B 
by a number of methods; we calculated the ratio of thermal to magnetic pressures, 
we calculated the critical field necessary to completely support the cloud against 
self-gravity and compared it to the observed field, and we calculated the ratio of 
mass to magnetic flux in terms of the critical value of this parameter. All these 
methods indicated that the magnetic field in S88B is dynamically significant, 
and should provide an important source of support against gravity. Moreover, 
the magnetic energy density is in approximate equipartition with the turbulent 
energy density, again pointing to the importance of the magnetic field in this 
region. 
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INTRODUCTION 



Massive stars affect the structure and dynamics of galaxies, the planet formation process, 
and even the formation of low mass stars via outflows, winds, and supernova explosions 



that drive mixing and turbulence in the interstellar medium (e.g., iZinnecker Sz Yorkei 120071 . 
and references therein). Yet, our understanding of their formation and evolution remains 
unclear. Moreover, the importance of magnetic fiel ds in the formation of stars, both low 



and high mass , has long bee n acknowledged (e.g., iMouschoviasI Il987t iHeiles et al.l 11993 



Crutcher 1999: McKee 



19991 ). Indeed, ongoing theoretical work continues to probe the 



effects of magne t ic fields on the star formation process (e . g.. IZhao et al. 



: Boss 


2007; 


Shu et al. 


2006: 



2011 



Mellon fc Li 



2005). However, 



observational data on magnetic fields is still scarce and theoretical work is hampered by the 
lack of information on magnetic field strengths in different density regimes. 

With these considerations in view, we have observed the high mass star forming region 
S88B with the Karl G. Jansky Very Large Array (VLA) for the Zeeman effect in thermal 
absorption lines of OH at 1665 and 1667 MHz. Observations of the Zeeman effect in absorp- 
tion lines with interferometers like the VLA provide an excellent method of mapping the 
magnetic fiel d in regions that are along the line of sight tow ard strong background continuum 



sources (e.g., iBrogan fc Trolandll200ll : ISarma et al.ll2000l ). In the next subsection, we will 
present some details of the source S88B. In § [2J we present details of the observations and 
reduction of the data. The results are presented in § [31 and discussed in § HI In § we state 
our conclusions. 



1.1. S88B 



S88B is a high mass st ar forming region (n ear I = 61?48,6 = 0?09) with a bolometric 
luminosity of 9.0 x 10 4 L© (IMueller et all 12002 1. It has been observed at a wide variety of 
wavelengths. There is some confusing nomenclature in the literature, and so we begin by 
discussing the morphology of its surroun dings. S88 is itself a diffuse extended nebula, and 



is likely a very evolved H II region (e.g., iFelli fc Hartenlll98ll ). About 14' to its southeast 



lie two bright knots of o ptical nebulosity: S88A and S88B (e.g., see the red Palomar Sky 
Survey print in Fig. 1 of lEvans et al.lll98ll ). No radio continuum emission is observed from 
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the western source, S88A. Strong radio emission is observed from an H II region centered 
slightly to the east of the optical ne b ulosit y S88B, and this is the target source for our 
Zeeman observations. iFelli fc Hartenl ( 119811 ) first proposed that the radio source in S88B 
was composed of two components, based on their 5 GHz continuum observations. The more 
extended component to the west (e.g., see § I3.1[ and Fig. 1 in this paper) that intersects 
the optical nebulosity is S88B-1, whereas the compact component to the ea s t is S 88B-2. In 
their high resolution 6 cm VLA image (HPBW 0.4"), IWood fc Churchwelll ( 119891 ) resolved 
out many details of these two sources — what they call G61.48+0.09A is referenced as 
S88B-2 in this paper (and most of the literature), whereas their sourc e G61.48+0.09B is a 



group of low intensity sources to the southwest of S88B- 1. Following lEvans et al.l ( 119811 ). 
we adopt the distance to S88B of 2.0 kpc. More recently, ISewilo et al.l ( 120041 ) have claimed 
a distance of 4.1 ± 2.4 kpc to S88B on the basis of their HllOa recombination line and 
H 2 CO abs orption obse r vation s; we note that our adopted value is within their error limits. 
Moreover, iGoetz et al.l ( 120031 ) have estimated on the basis of foreground star counts that the 
distance to S88B cannot be substantially greater than 2.0 kpc. 



2. OBSERVATIONS & DATA REDUCTION 

The observations were carried out with the Very Large Array (VLA) of the NRAO0 
in 2003 in the B-configuration, and combined with C-configuration data observed in 1997. 
Important parameters of the observations are listed in Table [H Both circular polarizations 
and both main lines (1665 and 1667 MHz) of OH were observed simultaneously in absorption. 
In order to mitigate instrumental effects, a front-end transfer switch was used to alternate the 
sense of circular polarization passing through each telescope's IF system every ~10 minutes. 

The editing, calibration, Fourier transformation, deconvolution, and calculation of opti- 
cal depths for the OH data were carried out using the Astronomical Image Processing System 
(AIPS) of the NRAO. Following standard procedure for Zeeman effect observations, band- 
pass correction was applied only to the Stokes / data set, since bandpass effects subtract out 
to first order in the Stokes V data, and the bandpass correction process adds noise to the 
V profiles. The B and C configuration data were combined in the uv-plane using the AIPS 
task DBCON . The rms in the continuum image of the combined B and C configuration 
data is 0.8 mJy beam -1 , whereas the rms in a channel of the line data is 1.7 mJy beam -1 . 
Magnetic fields were determined as described in § 13.31 by using the Multichannel Image Re- 



1 The National Radio Astronomy Observatory (NRAO) is a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, Inc. 
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construction, Image Analysis and Display (MIRIAD) software processing package, formerly 
of the Berkeley-Illinois-Maryland Array (BIMA) and now incorporated into the Combined 
Array for Research in Millimeter- Wave Astronomy (CARMA) telescope. 



RESULTS 



3.1. The 18 cm Continuum 



Figure [T] shows our highest resolution 18 cm continuum image of S88B made with 
uniform weighting (beam HPBW 3'.'1 x 3'.'0). The extended western component S88B-1 an d 
the compact component on the eastern side S88B-2, first identified by iFelli fc Hartenl (|198ll ). 
are clearly visible in this figure. We determined the integrated flux for S88B-1 to be 4.27±0.23 
Jy, and the integrated flux for S88B-2 to be 0. 79 ± 0.32 Jy. In a 5 GH z continuum survey 
with the NRAO 140-ft telescope (HPBW 6f5). iReifenstein et all Jl97oh obtained 6.1 Jy for 
the source 61.5+0.1 which is a combination of S8 8B-1 and S88B-2. In a 12-hr synthesis map 
with Westerbork (WSRT) at 5 GHz (HPBW 6" 6l.lFelli fc Hartenl dl98lh obtained a total flux 
density of 6.1 ± 0.2 Jy for S88B. At 4.9 GHz, iGarav et all fjl993h measured 4.4 Jy toward 
S88B-1 and 0.8 Jy toward S88B-2 in a combined C and D configuration observation with the 
VLA (HPBW ~ 4"). The integrated flux densities measured with these shorter wavelength 
data are in good agreement with the 18 cm VLA results, suggesting that the emission is 
predominantly either near the turnover point or optically thin (a = —0.1, S u oc u a ). Figure 



[2] shows the 350 /j,m dust continuum image of iMueller et al.l (120021 ) superposed on our 18 cm 
continuum image, along with other pertinent information — these are discussed in § 14.11 



3.2. OH Absorption 



The OH absorption profiles toward S88B are complicated, with several blended compo- 
nents toward a number of positions. In order to characterize the velocity structure, optical 
depth profiles for t he OH lines were calculated with the AIPS task COMB using standard 
procedure (e.g., see lRoberts. Crutcher. fc Trolandlll995l ). The OH column density was then 
determined using the relation 



Nqk/T cx = C r v dv 



cm 



K" 



T 



where T ex is the excitation temperature of OH and the constant C = 4.1063 x 10 14 and 



2.278 5 x 10 14 cm 2 (km s 1 ) 1 for the OH 1665 and 1667 MHz lines respectively (I Crutcher 

mi)- 
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The optical depth profiles reveal complex kinematics with a number of heavily blended 
components present toward much of the region. In order to highlight the prominent ab- 
sorption components, contour plots of the OH optical depth in several velocity channels are 
shown in Figure El superposed on a grayscale image of the 18 cm continuum. Every third 
velocity channel is displayed, and the plots are shown over the velocity range 17.4 km s~ x to 

24.0 km s" 1 . Figure H shows plots of iV(OH)/T cx toward S88B at 1667 MHz. In principle, 
variations in iV(OH)/T cx may arise due to variations in the OH column density, or due to 
variations in T ex , or both. However, in the results and discussion, we quote only variations 
in iV(OH)/Te X , since the excitation temperature cannot be measured based on absorption 
studies alone. Further discussion about possible variations in T cx appears in § 14.11 

From Figures [3] and HI we identified four prominent absorption enhancements, and these 
are marked in Figure H] as OH-a, OH-b, OH-c, and OH-d respectively. The absorption 
enhancement OH-a lies to the southwest of the continuum peak of S88B-1; it is seen in the 

19.1 km s _1 panel in Figure [31 with a peak r ~ 0.5. OH-a is also seen in the 19.9 km 
s _1 panel of Figures [3] with a higher peak r ~ 1. In this 19.9 km s" 1 panel, we also see 
the enhancements OH-b and OH-c, the former to the northwest of the continuum peak of 
S88B-1, and the latter about midway between S88B-1 and S88B-2. Finally, OH-d is seen in 
the 20.7, 21.5, and 22.4 km s _1 panels of Figure [31 with high peak r ~ 2 in two of these 
panels; it lies northward of the S88B-2 continuum peak; OH-d is also visible at the highest 
velocities at which OH absorption is observed (~ 23 km s _1 ). The high r's in OH-d account 
for the peak in iV(OH)/T cx to the north of the continuum peak of S88B-2. 



3.3. Magnetic Field Strengths 



The magnetic field strengths in S88B were determined using the Zeeman effect, by fitting 
a numerical frequency derivative of the Stokes I spectrum to the Stokes V spectrum for each 
pixel in the absorption line cube. The Stokes parameters have the standard definition of I = 
(RCP + LCP), and V = (RCP - LCP); RCP is right- and LCP is left-circular polarization; 
RCP has the standard radio definition of clockwise rotation of the electric vector when viewed 
along the direction of wave propagation. Since the observed V spectrum may also contain a 
scaled rep lica of the I spectrum itself, the Stokes V spectral profile can be described by the 
equation (ITroland fc Heileslll982l ; ISault et al.lll990l ): 



V =al + 



b dl 
2 d7 



(2) 



In this equation, the fit parameter a represents the scaled-down replica of the I profile in the 
V profile; a is usually the result of small calibration errors in RCP vs. LCP; for all results 
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reported here, a < 10 -3 . For thermally excited lines such as the OH 1665 and 1667 MHz 
lines reported in this paper, b = zB cos 9, where B is the magnetic field, 9 is the angle of 
the magnetic field to the line of sight, and z is the Zeeman splitting factor which is equal 
to z = 3.27 Hz nG- 1 for the 1665 M Hz line of OH, and 1.96 Hz /xG" 1 for the 1667 MHz 



line of OH ( iKazes fc Crutcherl Il986l ). We determined magnetic field strengths in S88B by 
using a least- squares method to fit equation (jSJ) to the Stokes V spectra to solve for a 
and b. The results of the fits give the line-of-sight component of the magnetic field, Bi os 
= B cos 9. It is usual practice to consider the results to be significant only if the derived 
value of -Bi os is greater than the 3a level. For Zeeman effect observations in OH absorption 
lines, we usually impose a stronger condition — we consider the detections to be significant 
only if both 1665 and 1667 MHz B\ os values are greater than the 3a level. In the current 
observations, however, this stronger criterion for significance would leave out some important 
information on the field. Toward S88B-2, the 1665 MHz profiles show a ~200 fiG field with 
about 5a level of significance, whereas the 1667 MHz profiles show a field below the 3a level. 
Therefore, for S88B-2 only, we have relaxed the criterion above to consider the measurement 
to be significant if it displays a 3a detection only in the 1665 MHz line. We believe this 
is reasonable, since the measured value of B\ os toward this source does not appear to be 
wildly discrepant from the measured value toward S88B-1. The non-detection at 1667 MHz 
toward S88B-2 is not surprising. The 1667 MHz Stokes V signal is weaker than the 1665 
MHz Stokes V signal because the Zeeman splitting factor for OH 1667 MHz lines (z = 1.96 
Hz /iG _1 ) is less than that for OH 1665 MHz lines (z = 3.27 Hz /iG -1 ). This inequality in 
the Zeeman factors would not impact the result as much if the OH gas was thermalized and 
homogenous, since the intensity ratios would then be r 67 /r 65 = 9/5. However, as described 
below (§ I4.ip . tqj/tq 5 = 1.0 — 1.3 toward S88B, which makes the sensitivity of the 1667 MHz 
line to the Zeeman effect less than that of the 1665 MHz line. 

The resulting B\ os map (taken from the 1665 MHz data) is shown in Figure [5J In the 
area enclosed by the thick-lined ellipse in this figure, the detected field is above the 3a level 
only in the 1665 MHz line, whereas in the rest of the displayed B\ os map, the field is above the 
3a level in both 1665 and 1667 MHz lines, with peak signal-to-noise ratios of ~10. Examples 
of the Stokes / and V profiles together with the derivative of / scaled by the fitted value of 
i?i os ar e shown in Figure M and Figure [7J the locations toward which these profiles have been 
plotted are marked in Figure [5] by a "+" and a "x" respectively. From Figure |6] and Figure 
[7J it is evident that the magnetic field is detected in the component with ^lsr ~ 21 km 
s _1 ; no fields were detected in any of the other velocity components. Toward S88B-1, we see 
from Figure [5] that -Bi os increases from about 90 fiG in the southwest to about 400 /iG in the 
northeast; this is the case in both 1665 and 1667 MHz lines. Meanwhile, toward S88B-2, the 
1665 MHz line reveals values of B\ os as high as ~270 /xG. 
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DISCUSSION 



4.1. Morphology and Kinematics 



Early optical, radio, and infrared observations of S88B by iPipher et al.l (119771 ) revealed 
that the peak of the near infrared (NIR) position is eastward of the optical nebula, whereas 
the radio continuum peak is even farther to the east than the NIR peak. This eastward 
progression is consistent with our observation of the location of the absorbing OH gas. It is 
clear from Figure [3] that almost all of the OH absorbing gas is located in the eastern half of 
S88B (i.e., all of S88B-2, and the eastern part of the extended source S88B-1). We did not 
detect any OH absorption in the western parts of S88B-1, where it is coi ncident with t he Ho ? 
nebulosity. The molecular cloud associated with S88B was mapped by Evans et al.l (119811 ) 
in several CO transitions at low resolu tion (HPBW 1/-2'), and found to extend over 10' (6.5 
pc) with a total mass of 5 x 10 3 M . lEvans et al.l (Il98ll ) speculated that the H II regions 
S88B-1 and S88B-2 are a blister formation in the molecular cloud which we are observing 
from the side. The velocities of the absorbing OH gas matc h up well with m olecular line 
velocities; e.g., the CO lines are si ngly peaked at 21 km s _1 ( lEvans et al.l 1 198 ll ). Moreover, 
recombination line observations by lGaray et al.l ( 119941 ) with the VLA (HPBW ~ 2") revealed 
that the center velocity of the single-peaked H92a line emission toward S88B-1 increases from 
22 km s -1 at the eastern edge of S88B-1 to 32 km s" 1 at the western edge. The 22 km s" 1 
velocity in the H92a recombination line at the eastern edge of S88B-1 matches well with the 
OH absorption data — e.g., see the panels between 21.5 km s _1 and 22.4 km s _1 in Figure EJ 
On the other hand, the recombination line velocit y at the western edge matches well with 
the optical Ha velocities which exceed 28 km s _1 (jDeharveng fc Maucheratlll978l ). In other 
words, we are observi ng an edge-on transition from the ionized gas to the molecular gas in 
S88B, similar to M17 farogan fc TrolandlboOll ). 



While the velocity structure of OH absorption observed by us toward S88B is compli- 
cated, the highest values for N(OH)/T ex are to the north of the continuum peak S88B-2. 
While this may, in principle, be due to variations in N(OH) or T cx or both, other tracers 
indicate that the column of absorbing gas is indeed greater in the e a stern parts of S88B (i.e., 
toward S88B-2). Toward the Ha peak in the west, iFelli fc Hartenl (Il981j) found a relatively 
low A v = 5 mag. Based on their observations of IR fine structure lines, iHerter et al.l (119821 ) 
measured A v = 26 =b 9 mag at the posi tion of the 2.2 NIR /im peak. This 2.2 fim peak 
was resolved by iDeharveng et al.l ( 120001 ) into a system of three stars (82, 83, 84) located 
on the eastern edge of the optical ne bulosity — their po sitions are marked in Figure [2] In 
contrast, toward S88B-2 in the east, iGoetz et al.l ( 120031 ) found A v ~ 70 mag by comparing 
their observed Bra flux with 5 GHz radio continuum data from the literature. S88B-2 also 



coincides with the 350 /im dust continuum peak from iMueller et al.l (120021 ). as shown in 
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Figure 2. There is also an extension in the 350 /zm dust continuum image to the northwest, 
which coincides with the location of S88B-1 (Fig. 2). In fact, due to the high obscuration 
toward the eastern parts of S88B, a complet e census of the exciting/ ionizing stars of S88B-1 
and S88B-2 has proved impossible so far. iDeharveng et al.l (120001 ) suggested that star 82 
could be the ionizing/exciting source of S88B-1, based on its high luminosity and central 
location (see Fig. |2J), but admitted the likelihood that other stars also contribute to the 
excitation/ionization in the region. Other possible ca ndidates for exciti ng sources are the 
stars LI (for S88B-2) and L2 (for S88B-1), detected by lPuga et al.l (120041 ) at a near-infrared 
wavelen gth of 3.5 p , m (L'- band) and marked in Figure [2J however, the issue remains open 
because iPuga et al.l (120041 ) could not determine the spectral type of LI or L2 as they did not 
detect them in their J, H, or K' near-infrared bands. 

If the absorbing OH gas is thermalized and homogenous (i.e., not clumpy within the 
observing beam), the ratio of the optical depths of the 1665 MHz and 1667 MHz lines should 
be tqi/tq^ = 9/5 = 1.8. However, the observed profiles toward S88B do not show this ratio. 
Instead, observed ratios tqj/t^ are in the range 1.0-1.3 over most parts of the source. While 
it is likely that the excitation temperatures for the two lines at 1665 MHz and 1667 MHz are 
not the same, it is unlikely that this difference would cancel out the expected ratio of optical 
depths over the whole source. More likely, the OH gas is clumped on size scales smaller than 
the beam. Such dumpiness has been found also in molec ular emission lines. The C O J=2— 1 
and 13 CO J=2-l and J=l-0 data (HPBW 22"-33") of Iwhite fc Fridlundl ( jl992h reveal a 
horseshoe-like structure surrounding the optical nebulosity, which they ascribe to material 
excited along the pe r iphery of an outflow cavity that encloses the optical nebula. Indeed, 
Phillips &: Mampasol ( 119911 ) have reported a low collimation bipolar outflow extending ove r 
3' from their CO observations; t his CO outflow was also o bserved by lRidge &: Moord ( 120011 ). 
Both the CO and 13 CO images in I White &: Fridlundl (119921 ) show highly fragmented structure, 
in agreement with the clumpy structure suggested by the OH absorption observations. More 
evidence of clumping c omes from observati ons of NH 3 lines at high angular resolution (HPBW 
4") with the VLA by I Gomez et al.l (119951 ). They found compact ammonia structures with 
sizes of 0.2 pc toward both S88B-1 and S88B-2. They suggest that these ammonia clumps 
correspond to compact molecular structures embedded within the larger molecular cloud, 
and that these clumps have been heated by the radiation fro m the star( s ) that ionize(s) the 
associated H II regions. Finally, recent C 18 observations by lSaito et al.l ( 120071 ) with a beam 
HPBW of 15" provide additional evidence of clumpy structure (also see § 14.2. II below); their 
clump B is marked in Figure [2j 

The OH excitation temperature T ex (which will be required in the calculations below) 
cannot be measured based on absorption studies alone, but reasonable limits for it can 
be set from our observations, and a likely value estimated based on other tracers. The 
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minimum continuum S u against which the line is detected is ~ 10 mJy beam -1 — we can 
use this to set an upper limit on T cx . Using S u — > : 1 mJy beam -1 —)- 17 K, we get 
an upper limit of 170 K for T ex . It is also unlikely that T cx is less than the typical dark 
cloud value of 10 K. Estimates of the kinetic temperature from other molecular tracers 
are also relevant if we assume LTE. From observations of the (2,2) and (3,3) lines of NH3, 
Gomez et al.l (119951 ) found Tkm = 80 K toward S88B-1 (their western N H3 clump) and 70 



K tow ard S88B-2 (their easter n NH3 clump 

( 



. Bas ed on CO observations, IWhite fc Fridlund 
(jl992l ) estimated T kiT1 ~ 60 K. IChini et al.l (Il986h f o und a dust temperature for S88B equal 
to Tdust = 42 K. iMozurkewich. Schwartz, fc Smith! (119861 ) fitted a 2-component dust model 
for S88B, and obtained T^ust = 40 K and 116 K. The higher may be du e to compact 
embe dded sou rces or it may reflect an enhanc ed population of small grains (jWu fc Evans 



19891 ). Indeed, iJourdain de Muizon et al.l (119901) have reported observations of PAH emission 
features toward S88B. iMueller et al.l (120021 ) found T dust = 75 K from their 350 /im dust 
continuum observations. We adopt T cx = 60 K for our calculations, but note that T ex likely 
varies over the source; still, it must lie between 10-170 K. 



4.2. Magnetic Fields 



Our observations of the Zeeman effect in OH absorption lines at 1665 and 1667 MHz 
have yielded B\ os values in the range of 90 to 400 /iG. These values are typical of magnetic 
fields in molecular clouds detected via the Zee man effect in the density regimes accessible 
via H I and OH thermal lines, as tabulated in ICrutcherl (11999). Note, however, th at our 
detected fields are higher than the single dish measurements by ICrutcher et al.l (119871 ) , who 
found B\ os = 69 ± 5 /iG toward S88B with the Nancay single dish. This is expected, since 
the low angular resolution (HPBW 3'.5 x 19') of Nancay will cause fields to be averaged over 
a large area, leading to much lower B\ os values. 



4-2.1. Magnetic Field Energetics 

The magnetic field in a cloud may consist of a static component B s , and a time- 
dependent or wave component B v . The static component connects the cloud to the external 
medium and determines the total magnetic flux throughout the cloud, but it can only pro- 
vide support to the cloud perpendicular to the field lines. The wave component is associated 
with MHD waves in the cloud, and it can provide three dimensional support to the cloud. In 
principle, the observed B\ os may consist of contributions from each of these two components. 
In reality, however, the Zeeman effect likely traces the strength of the static component of 
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the magnetic field only. This is because spatial averaging over the beam and along the line 



of sig ht tends to cancel the contribution from the wave component (e.g., iBrogan fc Troland 
200ll . and references therein). A principal goal of Zeeman effect measurements, therefore, 
is to estimate the importance of the magneti c field to t he dy namics and evolution of star 
forming regions like S88B. As summarized in ICrutcherl (119991 ). such estimates also require 
other physical parameters such as internal velocity dispersion, hydrogen column density, and 
radius, to be known. 

The relevant physical parameters mentioned above for estimating the importance of 
the magnetic field are available either from the present observations of S88B or from the 
literature. First, we must decide on an appropriate value for the radius. The highest res- 
olution molecular line emission observations available to date are the C 18 observations of 
Saito et al.l (120071 ) referred to above (with beam HPBW 15"). Clumps B, C, and D in their 
observations all have radius ~ 0.2 pc. Therefore, we adopt r = 0.2 pc for the radius in 
all the following calculations. Next, the hydrogen column density can be estimated from 
the OH line data, subject to certain assumptions. We use an average value of Non/T ex 
= 6.0 x 10 14 cm -2 K" 1 (see Fig. @J. For the excitation temperature, we use T cx = 60 
K, based on the discussio n in § I4TT1 abov e . Th en, to get Nn 2 , we use the conversion ratio 
NkJNqu = 2 x 10 6 from iRoberts et al.l dl995h for SI 06 This makes th e OH abundance 
a factor of 10 greater than that in dark clouds found by ICrutcherl ( 119791 ). Note, however, 
that if we consider the standard conversion ratio N-& 2 /A w ~ 10 21 , and use A v = 70 mag 
for S88B (see § 14.11 above), we obtain Nh 2 /Nqh = 1.9 x 10 6 , which justifies our use of the 



Roberts et al.l (119951) value rathe r than the dark cloud value; similar results were found for 
M17 by lBrogan fc Trolandl (120011 ) . Using this conversion ratio, we get Nn 2 = 7.2 x 10 22 cm -2 . 
This gives a hydrogen density n^H^) = 1.7 x 10 5 cm" 3 , w hich compares well with the value 



of n(H 2 ) = 1.51 x 10 5 cm 3 found by ISaito et al.l (120071 ) for the clump designated as B in 



their C 18 Q observations (se e Fig. El) . It is worth noting here that at lower resolution (HPBW 
30"), iPhillips fc Mampasol ( 119911 ) determined a mean n(H 2 ) ~ 6 x 10 3 cm" 3 . Together, 



this is further proof of the clumping in S88B discussed in 14.11 above, since the larger beam 
filling factor in a higher angular resolution observation of a clumpy region will naturally 
give a larger value for the density. Note, howeve r, that the density c annot be substantially 
greater than ~10 5 c m" 3 , a conclusion reached by lPhillips et al.l (119881 ) based on the absence 
of H2CO emission (jEvans et al.lll98ll ) toward S88B. Other parameters of interest are the 
velocity dispersion (which can be obtained from the observed line width) and the magnetic 
field; the adopted values for the average observed line width of OH lines and the averag e 
B\ os in S88B are listed in Table [2J For the calculations below, we will follow ICrutcherl ( 119991 ). 
and use total (static) magnetic field strength equal to 2 times Bi os . 



We can now proceed to assess the importance of the magnetic field in S88B, using a 
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number of approaches. First, we con sider the parameter /X,, which is the ratio of thermal 



to magnetic pressures (ICrutcherl 119991 ); it is a crucial parameter in any theory or simulation 
of the structure and evolution of molecular clouds that incorporates magnetic fields. If 
p < 1, then magnetic fields are important; the lower the value of p , the more dominant 
is the magnetic field. Based on the adopted values in Table [2j we find (3 P = 0.14 in S88B, 
indicating the importance of the magnetic field in this region. Our value for (3 P in S88B is 
similar to that in high m ass star forming regions like DR 21 OH 1, where (3 P = 0.21, and W3 



main), where (3 P = 0.13 (ICrutcherlll999l . and references therein). Another way to assess the 



importance of the magnetic field in S88 is to use the relation 

B s , crit = 5 x 10^ 21 N p /iG, (3) 

where N p is the average proton column density in the cloud. Equation[3]is a result of equating 
the static magnetic energy of the cloud to its gravitational energy; Bs, crit is then the average 
static magnetic field in the cloud that would completely support it against self-gravity. If 
the actual static magnetic field in the cloud B s > B s , CT it, then the cloud can be completely 
supported by the magnetic field, and is said to be magnetically subcritical. Further evolution 
of the cloud perpendicular to the field lines will occur primarily via ambipolar diffusion. 
The field can be judged to be dynamically important to the region even if it is less than 
Bs, crit; but is comparable to it. For Nn 2 = 7.2 x 10 22 cm -2 from Table [2] above, we get 
Bs, crit = 700/iG. Our average total (static) magnetic field strength (= 2 times B\ os in 
Table [2]) is 400 /iG; that is, the observed magnetic field strength is within a factor of ~2 
less than the critical field. Therefore, the magnetic field should be dynamically significant, 
providing an important source of support against self gravity. Yet another parameter to 
judge the extent to which a static magnetic field can support a cloud agai nst gravit a tiona l 



collapse is the observed mass to magnetic flux ratio [M/$ B ], defined by ICrutcherl (119991 ) 
in units of the critical value. This is equivalent to the consideration of Bs. C rit, an d may 
be redundant here, but we state it nevertheless to express our results in this alternative 
language. If the ratio of observed [M/$b] to the critical value is greater than 1, then the 
region is magnetically supercritical, meaning that the magnetic field cannot by itself provide 
support against gravitational collapse. For S88B, we find the observed [ M/$r1 in units o f 



the critical value to be 1.8. This is consistent with observations to date; ICrutcherl (119991 ) 
found the observed [M/$b] in terms of the critical value to have an average value equal to 2 
in a summary of available Zeeman measurements in molecular clouds. The calculated value 
for S88B indicates that whereas the static magnetic field by itself is not enough to support 
the molecular cloud against gravitational collapse, it is high enough to provide a significant 
source of support. Perhaps more important, however, in a highly supersonic medium like 
the S88B molecular cloud, is that for the magnetic field to be important, its energy density 
must be comparable to, or greater than, the turbulent energy density. The supersonic nature 
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of the medium is clear from the value of the sonic Mach number m s given in Table 3; from 
it, we have that the ratio of the one-dimensional velocity dispersion (a) to the isothermal 
sound speed (c s ) in S88B is u/c s = 2.3. This is a common result in molecular clouds; in 



his compilation of results for 27 molecular clouds, ICrutcherl (119991 ) found that motions are 
supersonic by about a factor of 5. Table 3 also lists the Alfvenic Mach number (m^) to be 
equal to 1. This implies that the magnetic energy density is comparable to the turbulent 
energy density, that is, these two energies are in approximate equipartition. 

Finally, we can use the measured B\ os to assess the relative importance o f the gravita 



tional, kinetic, and magnetic energy terms W, 7~, and A4, in the virial equation (IMcKee et al. 



19931) . To calculate their values, we use the expressions for W, T, and Jv[ given in ICrutcher 



1999|); the values are given in Table El We find that M « (1/2)T « (1/4) W, again in- 
dicating that magnetic fields are a significant, but not the sole, means of support against 
gravitational collapse. Finally, we find that 2T + M. = 1.3W; however, the surface pressure 
term in the virial equation, which is not included here, will act in the same sense as gravity; 
therefore, S88B is in approximate virial equilibrium. 



5. CONCLUSIONS 

We have observed OH thermal lines in absorption at 1665 and 1667 MHz toward the 
star forming region S88B with the aim of mapping magnetic fields via the Zeeman effect. 
The OH absorption line profiles toward this source are complicated, and contain several 
blended components toward a number of positions. The OH absorbing gas is located in 
the eastern parts of S88B, toward S88B-2 and the eastern section of S88B-1. The ratio of 
optical depths r 67 /r 65 = 1.0 — 1.3 over most parts of the source, instead of the value of 1.8 
for thermalized and homogenous OH gas. This is most likely due to clumping of OH gas, 
a conclusion supported by observations of CO, 13 CO, NH3 and C 18 emission lines from 
the literature. From the N(OH)/T cx plots, we identified 4 prominent OH clumps, which we 
have designated as OH-a, OH-b, OH-c, and OH-d. The highest values for N(OH)/T cx are 
to the north of the continuum peak S88B-2. Although this could, in principle, be due to 
an enhancement in N(OH) or T cx (or both), it is likely that north of S88B-2 is truly the 
location of the largest column of absorbing OH gas, as supported by independently measured 
extinction values toward S88B from the literature. Based on the observed quantities and 
adopted values of parameters from the literature, we obtain a gas density n(H 2 ) = 1.7 x 10 5 



cm 3 which compares well with the value of n{H-i) = 1.5 x 10 5 cm 3 found by ISaito et al. 
(120071 ) for a C 18 clump coincident with the S88B-1 and S88B-2 sources. 



In brief, S88B appears to have an overall geometry similar to Ml 7, in that there is 
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an edge-on progression from ionized to molecular gas going eastward. Our magnetic field 
map appears to mirror this transition, at least over the extended S88B-1 region, where we 
detected line-of-sight magnetic fields ranging from 90-400 fiG, with the lowest values to the 
southwest of the S88B-1 continuum peak, and the highest values to its northeast. Toward 
S88B-2, the fields are significant at the 3-a level in the 1665 MHz line only, and are detected 
over a small region but if we had the sensitivity to keep going eastward, the field would likely 
follow the same increasing trend. 

The detected fields allow us to assess the importance of the magnetic field in S88B, using 
a number of methods. We find that the ratio of thermal to magnetic pressures, P = 0.14; 
a value of /3 P < 1 means that magnetic fields are important; the lower the value of /5 P , 
the more dominant the magnetic field. A second method to assess the importance of the 
magnetic field is to compare the observed magnetic field with the critical field that would be 
required to completely support the cloud against its self-gravity. We find that the observed 
magnetic field is within a factor of 2 less than the critical field. A third method to assess the 
importance of the magnetic field, equivalent to the second, is to find the ratio of observed 
mass to magnetic flux in terms of the critical value; for S88B, we find this to be 1.8. These 
considerations lead us to conclude that the magnetic field in S88B is dynamically significant, 
and provides an important source of support against gravity. Moreover, the magnetic energy 
density is comparable to the turbulent energy density, implying that these two energies are 
in approximate equipartition, and pointing again to the importance of the magnetic field in 
this region. 
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Table 1. PARAMETERS FOR VLA OBSERVATIONS 

Parameter B- configuration C-configuration 

Date 2003 Nov 28/29, Dec 1/2 1997 Aug 15/16 

Configuration B C 

R.A. of field center (J2000) 19 ft 46 m 47 s .7 19 h 46 m 47 s .7 

Decl. of field center (J2000) 25°12'56".l 25°13'05".2 

Total bandwidth (MHz) 0.20 0.20 

No. of channels 128 128 

Channel Spacing (km s" 1 ) 0.28 0.28 

Approx. time on source (hr) 6.7 8.9 

Rest frequency (MHz) 1665.402 1665.402 

1667.359 1667.359 



Note. — l'=0.6 pc, for adopted distance to S88B = 2.0 kpc (see § II. ip 
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Table 2. Source Parameters 



Parameter Value 



Radius (r) 0.2 pc 

Column Density (7V H2 ) 7.2 x 10 22 cm" 2 

Particle Density (n p ) 1.7 x 10 5 cm -3 
Kinetic Temperature (Tk) 60 K 

Adopted Average FWHM Linewidth (Av) 2.5 km s _1 

Average Magnetic Field (B ios ) 200 yuG 
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Table 3. Source A : Derived Values and Virial Estimates 



Parameter Value 



m s 4.0 

rriA 1.0 

P p 0.14 

[M/<J> B ] obs/crit 1.8 

W 1.3 x 10 46 ergs 

T/W 0.54 

M/W 0.22 



Note. - m s is the sonic Mach number equal to \/3a/c s , where a is the line velocity 
dispersion and c s is the speed of sound; rriA is the Alfvenic Mach number equal to y/3a/v\, 
where is the Alfven velocity, (3 P is the ratio of thermal to magnetic pressures, M/$b is 
the mass to magnetic flux ratio, W is the virial gravitational energy, T is the virial kinetic 
energy, and M. is the yirial magnetic energy; expressions for calculating all these terms were 
taken from ICrutcherl (119991 ). 



19 h 46 m 50 s 49 s 48 s 47 s 46 s 

a (J2000) 

Fig. 1. — Contour image of the 18 cm continuum toward S88B. The beam is 3'.'1 x 3'.'0 
(P.A.=— 6?1), and the rms in the image is 0.8 mJy beam -1 . The contours are at 6, 12, 24, 
48, 96 mJy beam -1 . The extended western component S88B-1, and the compact component 
to the east S88B-2, are marked in the figure as Bl and B2 respectively. 
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Fig. 2. — Contour s and grayscale of the 350 fim dust emission toward S88B from 



Mueller et all (120021 ) superposed on the 18 cm continuum image. The HPBW of the 350 
/xm data is 14", and contours are at 1.4, 2 Jy beam -1 , and then in increments of 1 Jy 
beam -1 up to 10 Jy beam -1 . Our 18 cm continuum image has the same re solution and con- 



tours as in Figured] The positions of the NIR sources 8 2, 83, and 84 from iDeharveng et al. 



( 120001 ). the 3.5 um sources LI an d L2 from lPuga et all (120041 ). and clump B from the C 18 
observations of ISaito et all ( 120071 ) are also marked in the figure. 
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Fig. 3. — OH 1667 MHz optical depth channel images from 17.4 km s^ 1 to 24.0 km s" 1 
in ^lsr- The contours are from r = 0.25 to r = 2.0 in steps of 0.25. Every third velocity 
channel is displayed. The grayscale shows the 18 cm continuum at the same resolution as 
the OH optical depth image, which is 4'.'5 x 4'/3 (P.A.=14?1), as shown in the boxed inset in 
the top left panel. 
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a (J2000) 

Fig. 4. — Gray-scale image of NoH67/T ex toward S88B overlaid on the 18 cm continuum. 
The continuum image has the same resolution and contours as in Fig. [TJ The 1667 OH beam 
is shown in the lower right; its parameters are given in Fig. |3j 



-24- 



10 - 



13 00 - 



o 
o 

C\2 

•-3 

*0 



50 - 



40 - 



30 - 



25°12'20" h 



_i_ 




400 



300 



200 



100 



19 h 46 m 50 s 



49 s 



48 s 47 s 

a (J2000) 



46 s 



Fig. 5. — Image of the detected B\ os toward S88B, taken from the 1665 MHz data. The 
contours depict the 18 cm continuum, and are the same as in FigCD The thick-lined ellipse 
shows the region in which the detected B\ os is at the 3-er level in 1665 MHz only. The "+" 
sign indicates the position toward which the profiles are shown in Fig. [6] below, while the 
"x" indicates the position toward which the profiles are shown in Fig. [7j 
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V LSR ( km S ') V LSR ( km S ^ 

Fig. 6. — Stokes / (top; histogram) and V (bottom; histogram) profiles at 1665 MHz (left) 
and 1667 MHz (right) toward the position in S88B-1 marked in Fig. [5] by a "+" (a 2 ooo — 

19h 46 m 48 s 3 

^2000 = 25°12 / 48'. / 3). The curve superposed on V in the lower frame on the left 
shows the derivative of I at 1665 MHz scaled by i?i os = 159.4 ± 18.0 fiG, whereas the curve 
superposed on V in the lower frame on the right shows the derivative of / at 1667 MHz 
scaled by B los = 209.8 ± 29.5 fiG. 
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Fig. 7. — Stokes / {top; histogram) and V (bottom; histogram) profiles at 1665 MHz toward 
the position in S88B-2 marked in Fig. EJby a "x" (a 200 o = 19 h 46 m 49!0, (5 200 o = 25°12 / 44 / /7). 
The curve superposed on V in the lower frame shows the derivative of I at 1665 MHz scaled 
by -Bios — 272.6 ± 40.7 fiG; the fit for the field was made over a restricted range of velocity 
channels, and the channels over which the fit was made are indicated by the solid line. 



